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Abstract 

The exact matching condition is given for hadron matrix elements calculated 
in any two different schemes, in particular, in the lattice and dimensional reg- 
ularization, (modified) minimal subtraction MS schemes. The result provides 
insight into and permits to go beyond Lepage and Mackenzie's mean field 
theory of removing tadpole contributions in lattice operators. 
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Presently, lattice QCD provides the unique method with controlled approximation to 
compute hadron properties directly from the QCD lagrangian. In the last few years, a 
number of groups have calculated on the lattice an impressive list of hadron matrix ele- 
ments, ranging from the axial and scalar charges of the nucleon to lower-order moments 
of deep-inelastic structure functions [0-0] . Note, however, that most of the hadron matrix 
elements are not directly physical observables. In field theory, apart from the S-matrix, 
physical observables are related to symmetry generators of the lagrangian, such as the vec- 
tor and axial-vector currents or hadron masses. Nonetheless, hadron matrix elements are 
useful intermediate quantities to express physical observables. Being intermediate, they 
often depend on specific definitions in particular context. Or in field theory jargon, they 
are scheme-dependent. Since schemes are generally introduced to eliminate ultraviolet di- 
vergences in composite operators, the scheme dependence of a matrix element is in fact 
perturbative in asymptotically- free QCD. 

Understanding scheme dependence has important practical values. In calculating hadron 
matrix elements on a lattice, one is automatically limited to the lattice scheme. On the other 
hand, hadron matrix elements entering physical cross sections are often defined in connection 
with perturbation theory. The best scheme for doing perturbation theory is not the lattice 
QCD, because the lattice has complicated Feynman rules and accommodates only Euclidean 
Green's functions. The most popular scheme for perturbative calculations is the dimensional 
regularization introduced by t' Hooft and Veltman more than two decades ago, followed by 
the (modified) minimal subtraction (MS). 

A popular practice currently adopted in the literature for matching the matrix elements 
in the lattice and MS schemes goes like this [CM. Consider, for instance, a quark operator 
O. First, the one-loop matrix element of O in a single quark state \k) is calculated on the 
lattice, 

(^|0|A;) latt = A[l + ^( 7 lnay + C )] , (1) 



and in the MS scheme, 



(k\0\k) m =A^l + 9 -^( 1 \np 2 /^ + c')j , (2) 

where p 2 is an infrared cut-off. Second, the hadron matrix elements in the state \P) in the 
two schemes are assumed to have the relation, 

(P\0\P) Ws - (k\0\k)us~ + ^ ^ ^ () 

In the second equality, an expansion in g 2 is made. In order to cancel the infrared cutoff 
p 2 , one has to identity go(a) and g(fi). However, after the cancellation, one does not know 
which g should be used in the above equation. Due to the large tad-pole contributions, 
the difference between coupling constants in the two schemes is significant in the same 
momentum region. Furthermore, it is not clear how to generalize the above relation to 
multi-loops. 
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As indicated early, the difference between the matrix elements in any two schemes arises 
from the ultraviolet region. Thus there must exist an all-order perturbative relation between 
the matrix elements calculated in the two schemes. In this note, I shall derive such a relation 
using renormalization group arguments. 

Consider a physics observable M(Q 2 ) where Q 2 is a hard momentum scale of a physical 
process. M{Q 2 ) could be a moment of some deep-inelastic structure function. According to 
the operator-product expansion, we write, 

M(Q 2 ) = C(Q 2 / f , 2 ,a s ( f , 2 ))A( f , 2 ). (4) 

where C is the coefficient function, calculable in perturbative theory. C is clearly scheme- 
dependent, and let us assume a scheme has been chosen for its calculation. A is a soft 
hadron matrix element defined in the same scheme. The \i dependence of A satisfies the 
renormalization group equation, 

^ + 7 (^) = 0, (5) 

where 7 is the anomalous dimension of the operator in the scheme. The renormalization 
scale dependence must cancel between C and A, leaving M /i-independent. Therefore C 
satisfies the renormalization group equation, 

f-700C = 0. (6) 
The above equation can be integrated formally to yield, 

rs(Q) 7 



C(Q 2 /fi 2 , a s (fi 2 )) = C(l, a s (Q 2 )) exp ( - J 



Jg) , (7) 

where the QCD /^-function is also scheme- dependent. From the above formulas and the fact 
that M is scheme-independent, we derive a relation between the matrix elements in two 
different schemes, 

C(l, a s (Q 2 )) exp - / Idg A{y) = C(l, a s (Q 2 )) exp - / Idg A{fi) . (8) 

It seems that one has to know the coefficient functions in the two schemes before hand. 
However, this is unnecessary if one takes Q 2 — > 00 and uses asymptotic freedom. In the 
limit, the coupling constants in different schemes approach to zero at the same rate and 
C(1,0) becomes scheme- independent! Thus we have, 

This is the desired equation, which depends on the /9-functions and anomalous dimensions of 
the operator in two different schemes. The matrix elements defined in this equation satisfy 
the renormalization group equations in both schemes. 

To specialize to the lattice and MS schemes, let me make the following definitions, 
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where Os(l/e) and O^g are the bare and renormalized operators in the MS scheme, and 

r\ rrt \ r\ a dZ(a) 

Clatt(o) = Z(a)0 R l att ; 7latt = ^7 T ~, , (H) 

Z[a) da 

where Oi a tt(°) an d Or i att are the bare and renormalized operators in the lattice scheme. 
Define further a renormalization relation between the bare lattice and renormalized MS 
operators, 

ms(^) = Z(a s (a),a s (fi))Oi att (a) . (12) 
Then the renormalization constant is, 

Z(a 8 (a), a s (n)) = exp ( / ^-dg] exp ( - f -^-dg) . (13) 

yo Aatt / \ Jo %g J 

The corresponding relation between the matrix elements in the two schemes is Aj^({i) = 
Z{a s {a),a s {n))A utt {a). 

In practical applications, one uses a perturbative expansion of Eq. flT3|). To next to the 
leading order, 

7 *<£)' + .... 

" = -»( A( s )+A( s ),+ -) ' (14) 

Substituting the above into Eq. (|T3]) , I find, 



a. 



+ ... 



Z(a.(o),a a ( A *))=(^] 1 W/ - ^ I, (15) 

where I have used the fact that the first-two-loop (3 functions and the one-loop anomalous 
dimensions are scheme-independent. 71 is the two- loop anomalous dimension and is scheme- 
dependent. The constants c and d in Eqs. (1) and (2) do not enter in the above relation 
directly. However, they contribute through the two-loop anomalous dimensions. 

The first factor in Eq. ( f[5f) sums over all the leading logarithmic contributions. At the 
momentum scales associated with current lattice calculations, the strong coupling constant 
in MS is nearly a factor 2 larger than the bare lattice coupling because of the large tadpole 
contributions. Thus if 7o/2/5o is on the order of unity, the first factor in Eq. (15) is large. 
However, this is not the only source of large tadpoles contributions. The lattice operators, 
defined using the gauge link = exp^igoaA^), where is the gauge potential, also contain 
large tadpole effects, reflected by the large constant c in Eq. (1). This large constant induces 
a large two-loop anomalous dimension 71(a). As a consequence, the lattice perturbation 
series in the denominator of Eq. ([L5D converges sluggishly and must be improved with some 
resummation procedure. 
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In Ref. 0, Lepage and Mackenzie suggested to replace the bare lattice coupling a s (a) 
in a perturbative expansion by a renormalized coupling, such as the MS coupling, a s (fi). 
Such improvement procedure eliminates large coefficients in the expansion and speed up the 
convergence of the series. However, to completely cancel the tadpole contributions, one must 
multiply the series with an expansion of some power of a s (a)/a s (fi). The appropriate power 
is determined by the anomalous dimension of the operator, which in turn depends on a par- 
ticular construction of the operator on the lattice. The residual factor of a s (a) / a s (/j,) in Eq. 
(|15|) represents the intrinsic tadpole effect of the lattice operator, accounting for the major 
part of the renormalization factor. For the same continuum operator, one can construct 
many different versions of lattice operators, which are affected by tadpole contributions 
differently. 

In general, the lattice two-loop anomalous dimension 71(a) can give a good indication as 
to what power of a s (a) is needed to cancel the tadpole effects in the lattice series. Lepage and 
Mackenzie's mean field theory provides an empirical method to determine the residual power 
of a s (a)/a s (fi) in the renormalization factor. Consider a lattice operator constructed with 
k powers of gauge links. Then, according to Lepage and Mackenzie, the tadpole effects can 
be accounted for by multiplying the lattice operator with a factor of Uq k . Through studying 
the effect on the Wilson action, Lepage and Mackenzie determined u = (a s (a)/a s (/i)) 1//4 . 
Thus, the mean field theory predicts Z ~ (a s (a)/a; s (yu))~ A: / 4 apart from a well-behaved 
perturbation series. 

However, Eq. (pT5|) allows one to go beyond the mean field theory with the explicit depen- 
dence on perturbative corrections to high orders. With more powerful methods becoming 
available for computing lattice perturbation series, this exact formula for the renormaliza- 
tion factor is quite valuable. The formula also corrects the mean field prediction when the 
latter fails. An extreme case is the unit operator constructed as UU~ l on the lattice, for 
which a naive application of the mean field theory fails completely. 

To summarize, I have derived an exact relation between the matrix elements in lattice 
and MS schemes. The relation is used to discuss Lepage and Mackenzie's mean field theory 
for correcting the tadpole effects in lattice operators. The relation is useful for comparing the 
hadron matrix elements calculated in the lattice and measured, for instance, in deep-inelastic 
scattering. 
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